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Abstract 

We study the stationary states of a quantum mechanical system de- 
scribing an atom coupled to black-body radiation at positive tempera- 
ture. The stationary states of the non-interacting system are given by 
product states, where the particle is in a bound state corresponding 
to an eigenvalue of the particle Hamiltonian, and the field is in its 
equilibrium state. We show that if Fermi's Golden Rule predicts that 
a stationary state disintegrates after coupling to the radiation field 
then it is unstable, provided the coupling constant is sufficiently small 
(depending on the temperature). 

The result is proven by analyzing the spectrum of the thermal 
Hamiltonian (Liouvillian) of the system within the framework of W*- 
dynamical systems. A key element of our spectral analysis is the 
positive commutator method. 
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1 Introduction 



In this paper, we study a quantum mechanical model of an atom interacting 
with black-body radiation. The atom is described by an electron moving in 
a potential, e.g. the Coulomb potential of a static nucleus. 

Our goal is to show that when the atom is coupled to the quantized 
radiation field in the state of black-body radiation at sufficiently high tem- 
perature, it is ionized. We describe this ionization process by showing that 
stationary states of the system become unstable when the atom is coupled 
to the radiation field. Each bound state of the atom leads to a stationary 
state of the uncoupled system, where the field is in an equilibrium state. We 
consider a class of small interactions localized in space which couple finitely 
many bound states of the atom to the field. We show that the stationary 
states for the coupled system correspond to those eigenvalues of the atomic 
Hamiltonian which are not coupled to the field. In other words, all stationary 
states arising from atomic bound states which are coupled to the field by the 
interaction are unstable, provided the coupling constant is small enough. 

This instability is explained by the following mechanism. If the electron 
is in a bound state with energy E < 0, it will be hit, after some time, by a 
quantum (photon) of energy uo > —E, and hence will make a transition to a 
scattering state of energy E + uj. In other words, the atom is ionized. 

The average time, t E , it takes for an atom in a bound state of energy E to 
be ionized is given, to second order in the perturbation, by the Fermi Golden 
Rule. Heuristically, for an inverse temperature (5 of the field, it satisfies 

t E <xe® E \\p{E)\-\ (1.1) 

where |p(-E)| 2 is the probability for the electron to make a transition from 
the bound state to a scattering state by absorbing a photon of energy > E. 
The factor e® E \ in (1.1) can be explained by Planck's law, which says that 
the probability density for a photon to have energy oj is e pl_ x ■ At zero tem- 
perature, f3 = oo, one finds that t E = oo, and thermal ionization does not 
occur. For a given strength of the interaction, as measured by the size of a 
coupling constant A e 1, we are able to show that thermal ionization takes 
place, provided < |A| < ke~ 2 ^ E °^, where E < is the minimal energy of 
the electron in a bound state coupled to the radiation field, and k is some 
constant. This restriction is of technical nature; physically, thermal ioniza- 
tion is expected to be observed for arbitrarily small temperatures, provided 
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the coupling constant is small enough independently of (3. 

Next, we describe the system and our main results in some more detail. 
The atomic Hamiltonian is a Schrodinger operator H p = —A + v on the 
Hilbert space 7i p = L 2 (R 3 , d 3 x), where v belongs to a certain class of poten- 
tials including the Coulomb potential regularized at the origin. The operator 
Hp generates the Heisenberg dynamics 

a p t (A) = e itH *Ae- itH * 

on the von Neumann algebra 2l p = B{Ti p ) of bounded operators on H p . 

The field is conveniently described in terms of a C*-algebra 21/, which 
can be viewed as a time-averaged Weyl algebra. The dynamics is given by a 
♦automorphism group of 21/ describing free massless bosons. 

The combinded system is described in terms of the algebra 

21 = Stp <g> 21/, 

and the uncoupled dynamics is given by the automorphisms 

«t,o = of <S> a{. 

To define the coupled dynamics, we specify a (regularized) interaction term 
XV^ e \ whose form is motivated by standard models of atoms interacting 
with the radiation field. The regularization is introduced to guarantee that 
G 21, for all e ^ 0. The interacting dynamics, af[, is then defined as the 
♦automorphism group of 21 obtained by the Schwinger-Dyson series. 

At zero temperature, the dynamics of the model is generated by the formal 
Hamiltonian 

H = H P + H f + XV, 

where Hf = dr(|fc|) is the free- field Hamiltonian, i.e., the second quantized 
multiplication operator \k\, acting on bosonic Fock space JF(L 2 (1R 3 , d 3 k)), and 
the interaction term V is given by 

V = ^G a ®(a(g a ) + a*(g a )). 

a 

The sum is over a finite set, G a are bounded selfadjoint operators on B{7i p ), 
and the form factors g a are in L 2 (M 3 , d 3 k). 
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We introduce a reference state 

on 21, where u p is given by a (strictly positive) density matrix, and out is the 

(/3,a{)-KMS state of 21/, i.e., the state of black-body radiation at inverse 
temperature (5. We are interested in the time evolution of states on 21 which 
are close to (normal w.r.t.) u rei , i.e., which are represented by a density 
matrix on the GNS Hilbert space H of (21, u mi ). The GNS representation 
provides us with a representation map irp : 21 — > £>(7i) and a vector f2 rcf G 7i 
s.t. u/ cf (A) = (f2 ref , 7Tg(A)f2 ref ). There is a selfadjoint operator on 7i 
generating the coupled time evolution in the representation irp, 

M4l(A)) = e^ ) MA)e'^\ 
for all A e 21 and tel. We will show that 

s-lirne^' = e itLx 

exists, for all t, and defines a *automorphism group 

<r t ,\(A) = e itLx Ae- itL > 
of the von Neumann algebra 

Tip = 7173(21)" C B(H). 

The pair (971^, ct^a) is called a W^*-dynamical system. Our results concern the 
structure of the set of normal (a-weakly continuous) time-translation (<7t,\-) 
invariant states on QJTg. 

The general theory of von Neumann algebras shows that there is a one- 
to-one correspondence between normal a tt \-inva.Tiant states on OJtg and nor- 
malized vectors in the set 

VDkerLx, (1.2) 

where V is a certain cone in 7i, the so-called natural cone associated to 
(Tip, Q rci ), provided we choose the thermal Hamiltonian (Liouvillian) L\ in 
such a way that the unitary one-parameter group {e' ltLx | f G M} leaves V 
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invariant. Let M. be a labelling of the eigenvalues of H p , including multiplic- 
ities. An element m G M. is called a mode and the corresponding eigenvalue 
of H p is denoted by E(m). We will see that 

V n kerL = V H spanjy? m <g> y?„ ® | m,n E M,E(m) = -E(n) j, (1.3) 

where ip m is the eigenvector of H p corresponding to the mode m, and Q is 
the vector representative of LOp. Our main result, Theorem 2.3, shows that 
(1.2) is the subset of (1.3) with m, n ranging over those modes which are not 
coupled to the field. 

While this result holds for a specific class of potentials (see (2.2)), we 
prove in Theorem 2.2 a result which holds for a very general class of poten- 
tials: For any cr t)0 -invariant normal state a; and any cr^-invariant normal 
state oo x on 971^ we prove that \\lu° — lu x \\ > k > 0, provided A ^ is small 
enough, for a constant k independent of A. Here || • || denotes the norm on 
the space of linear functionals on Wlp. Theorem 2.2 is proven for bounded 
potentials v such that —A + v has only finitely many eigenvalues below the 
threshold of the continuous spectrum, all of which are coupled to the field. 
Alternatively, we could relax this finiteness condition but couple only finitely 
many modes to the field. 



2 Definition of the model and main results 

In Section 2.1 we introduce the model and show in which way it defines a 
W^*-dynamical system (%ftp,a tt \). Our main results are presented in Section 
2.2. 

2.1 Definition of the model 

Starting with an algebra 21 describing the joint system atom-field and a (reg- 
ularized) dynamics afl on it, we introduce a reference state u mi , describing a 
bound state of the atom and black-body radiation at inverse temperature f3. 
We then consider the induced (regularized) dynamics af\ on 7Tg(2l), where 
(H, Tip, £T ef ) denotes the GNS representation corresponding to (21, u Tcf ). As 
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e — > 0, (Tj^ tends to a *automorphism group, ct^a, of the von Neumann al- 
gebra Wlp, defined as the weak closure of 7Tg(2l) in B(H). We determine the 
generator, L\, of the unitary group, e ltLx , on H implementing a tt \] L\ is 
called a Liouvillian. We explain the relation between eigenvalues of L\ and 
invariant normal states on VJlp. 

2.1.1 Kinematical algebra 21, and regularized dynamics af\ 

We consider a system consisting of a quantum mechanical particle (an elec- 
tron in the potential of a static nucleus) interacting with a quantized field. 

Pure states of the particle system are described by unit vectors in the 
Hilbert space H p = L 2 (R 3 ,d 3 x), their dynamics is determined by the Schro- 
dinger equation with Hamiltonian 

H p = -A + v, (2.1) 

where the potential v is bounded and satisfies one of the following two con- 
ditions. 

• Condition Ca- The potential v is s.t. the spectrum of H p consists of a 
finite number d of eigenvalues (counting multiplicity) lying below the 
continuous spectrum which covers [0, oo). We set E := inf cr(H p ) < 0. 

• Condition Cb- The potential v is given by 

= -K^<1, (2.2) 

where p(|#|) is a smooth, non-negative function that has a zero of order 
1 + /i at the origin, and increases to a constant value p as \x\ — * oo, in 
such a way that v is smooth, and 

(x ■ V) j v (2.3) 

are bounded, for j — 0, . . . , 3. Notice that the eigenvalues of H p are all 
negative and can accumulate only at the threshold 0. 

Remark. In Condition Cb we admit potentials such that H p has infinitely 
many eigenvalues below zero, but couple only finitely many of them to the 
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field, as we explain below. 

The field is a scalar massless free bosonic field. (It would be more interest- 
ing, physically, to consider the quantized electromagnetic field. Our methods 
can be applied to the resulting model at the price of slightly more compli- 
cated notations.) The scalar free field is conveniently described in terms of 
a "time- averaged" Weyl algebra, 21/, which is the C*-algebra (of "observ- 
ables"), defined as follows. Let 22J be the Weyl algebra over the Hilbert 
space 

L 2 = L 2 (R 3 ,d 3 k)nL 2 (R 3 , |A;|-VA;), (2.4) 

i.e., 22J is the C*-algebra generated by Weyl operators, W(f), f G Lq, satis- 
fying the Weyl relations 

W(f)W(g) = e- ilm ^W(g)W(f). (2.5) 

The free field dynamics on 2U is given by the *automorphism group 

W(f) » afiWU)) = W(e^f), (2.6) 

where u(k) = \k\ is the energy of a single boson. For functions / G Lq, the 
expectation functional 

/ ^ e -i<^( 1+ ^i)^) (2.7) 

is well defined and determines a (f3, af°)-KMS state on 2U. It is well known 
that the *automorphism group af^ of 2U is not norm-continuous (i.e., R 3 
t i— > af°(W(f)) is not continuous in the norm of 2U). The time-averaged 
C*-algebra 21/ is generated by elements of the form 

a(h) = [ ds h(s)af(a), (2.8) 

where a G 22J and h : R — > C are functions whose Fourier transforms satisfy 
h G C£° (this is a convenient class of functions which allows us to define KMS 
states on 21/, see [FM]). The free field dynamics on 21/ is defined by 

a{(a(h)) = [ ds h(s - t)af(a) =: a(h t ). (2.9) 
Jr 
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It is a norm-continuous *-automorphism group on 21/. We refer to [FM] for 
more details on the construction and the properties of 21/. 

The joint system describing the particle and the field is described in terms 
of the C*-algebra 

21 = % <g> 21/, (2.10) 

where 2l p = B(7i p ) is the von Neumann algebra of all bounded opera- 
tors on the Hilbert space H p . The uncoupled dynamics is given by the 
♦automorphism group 

at,o = o? t <g> a{ (2.11) 

of 21, where c? t (•) = e ltHp . e ~ ltH p. In order to define the dynamics of the inter- 
acting system in a representation independent way (i.e., as a *automorphism 
group on 21), we need to introduce a regularized interaction term. For e ^ 0, 
this term is given by given by 

V # = ® ^- € { W ^9a)(h e ) - W(eg a )(h e y} G 21, (2.12) 

a 

where the sum is over finitely many indices a, with = e ^(W p ), 

<7 Q G Lq, for all a, and where fo e is an approximation of the Dirac distribution 
localized at zero. To be specific we can take h e (t) = \e~ t2 / t2 . The symbol 

G Qj # (and similarly V^) stands for either G a or G aj j, where J is some cutoff 
determining which modes of the particle are coupled to the field. In order 
to describe this more precisely, we introduce the following terminology. Let 
Ai be the index set of the discrete "modes" of H p , i.e., a labelling of the 
eigenvalues of H p including multiplicity. Given m G A4, E(m) denotes the 
corresponding eigenvalue of H p . An eigenvalue E of H p is simple if and only if 
there is a unique m G M. s.t. E = E{m). We denote the rank-one projection 
corresponding to the mode m G Ai by p m . 

Let Jd C M. be a set of finitely many discrete modes of H p and let J c 
be an open interval in the continuous spectrum M + of H p (we may also take 
a finite union of disjoint intervals), s.t. J c C [r,R], for some r,R satisfying 
< r < R < oo. The set 

J := J d U J c 
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determines the modes of the particle which are coupled to the field, according 
to the interaction 

G a ,j = (pj d + »(H p ))G a (p Jd + n(H p )), (2.13) 
where G a is a bounded, selfadjoint operator on Tt p , and 

pj d = E P"» ( 2 - 14 ) 

meJ d 

fi G C£°(J C ) is a smooth version of the indicator function with support in J c , 
and n(Hp) is defined via the Fourier transform 



Clearly, G a ^j tends to G a , in the strong sense as \x increases to the character- 
istic function of M + (i.e. J c f R + ) and increases to the set of all discrete 
modes of H p . Thus, V"j e ^ can be viewed as an approximation of = vj^ K . 

The interaction term (2.12) determines a *automorphism group af\ of 21, 
the coupled dynamics, via the norm-convergent Dyson series 

ft rtn-l 



rt rtn-i r 

ag(A) := a tfl (A) +J2(^) n dt ± --- dt n a tn , (V^), 

n>l J ° J ° L 

•••[a tl ,o(^ e) ),a ti0 (A)] •••]], (2.15) 

where A G 21, and A G M is the coupling constant. The multiple integral in 
(2.15) is understood in the product topology coming from the strong topology 
of 13 (Hp) and the norm topology of 21/. 

One may view af\ as a regularized dynamics, in the sense that it has a 
limit, as e — > 0, in suitably chosen representations of 21; (this is shown in 
[FM] and explained below). 

The functions g a G L\ are called form factors. Using polar coordinates 
in M 3 , we often write g a = g a (cu, E), where (cj, E) G IR+ x S 2 . 

We now specify two sets of assumptions on the interactions. 

Condition A. The potential v satisfies condition Ca, the interaction is 
given by V^ e \ and the following properties hold. 
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• Infrared and ultraviolet behaviour of the form factors. For any fixed E, 
g a (-, S) G C 4 (IR + ), and there are two constants < k±, k 2 < oo, s.t. if 
iq < ki, then 

\dlg a (uj, S)| < k 2 uj p - j , for some p > 2, (2.16) 

uniformly in a, j — 0, . . . , 4 and E G S* 2 . Similarly, there are two 
constant < Ki, K 2 < oo, s.t. if u > Ki, then 

\dlg a (uj, S)| < X 2 ^ _<?_j , for some q > 3. (2.17) 

• Relative bound on [G a ,H p ]. Define the commutator [G a ,H p ] in the 
weak sense on x by 

(V>, [G a , H p ]ip) = (G a ip, H p ip) - (Hptp, G a ip) . 

Then [G a , H p \ extends to a relatively (H p — _E + 1) ^-bounded operator, 
i.e. there is a k < oo s.t. for any *0 G Cq°, 

|| [G Q , ffJ^H < k \\(H P -E + lf^W , (2.18) 

where E = inf o~(H p ) < 0. 

• The Fermi Golden Rule Condition. We define a family of bounded 
operators on H p by F(u, E) = ^ Q (? a (u;, E)G Q and let, for arbitrary 
e > 0, 

T e ( W , £) = / dE E) — / c£ 2 F(^, S)*, (2.19) 

J S 2 {n p — h, — ujy + e z 

where is an eigenvalue of H p and p c is the projection onto the continu- 
ous subspace of H p . Let p(E) denote the projection onto the eigenspace 
corresponding to E. We assume that there is an e > 0, s.t. for 

< e < e , 



u 2 



duJ ~^—{P( E ) T ^ E )P( E ) ^ 7e p(E), (2.20) 
E e pw — 1 



for any E G a p (H p ), where 7^ is a strictly positive constant. We set 

7 := min 7^ > 0. (2.21) 

Ee<T P (H p ) 
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Remarks. 1) All requirements in Condition A are independent of the 
regularization of the interaction. 

2) For the physical model of an atom interacting with the radiation field, 
the value of the constant p in (2.16) is p = —1/2 (or p = 1/2 in the dipole 
approximation), see e.g. [BFS]. Although p > 2 is quite far from the phys- 
ical range, we do not attempt here to optimize condition (2.16). This will 
be the aim of subsequent work. Suffice it to note that the discrete values 
p — —1/2, 1/2 are also admissible in our analysis. 

3) The operator T e (u, E) is just a (non-negative) number if E is a simple 
eigenvalue. For e small, it represents the probability that the particle makes 
a transition from the bound state corresponding to the energy E into a scat- 
tering state with energy E + u > by absorbing a photon of energy uj. The 
probability density for a photon to have energy uj is given by Planck's law, 
i.e., by (e@ w — l) -1 . Hence 7 is a perturbative bound on the probability of 
an ionization process; it depends on the inverse temperature (5 as 7 ~ e^ E ° , 
where Eq < is the ground state energy of H p . More precisely, if we assume 
that, for < e < e , uj 1— > p{E)T e {uj , E)p(E) is continuous, and that there is 
a constant t > s.t. p(E)T € (uj, E)p(E) > t ■ p(E) at uj = —E, then one sees 
that fcfS < 1e < fce^, for some k which does not depend on (3. 

Condition B. The potential v satisfies condition Cb, the interaction is 
given by Vj € \ and the following properties hold. 

• The infra-red and ultra-violet behaviour of the form factors is as in 



where we set (x) = (x 2 + l) 1 / 2 , for x G M 3 . Notice that this is a 
condition on G a not depending on the regularization. 

• The Fermi Golden Rule Condition. For all eigenvalues E of H p s.t. 
E = E{m) for some m G Jd, let T e (u,E) be defined as in (2.19), with 
p c replaced by fi(H p ) 2 , and let pj d (E) = me.j d p m . There is an 



(2.16), (2.17). 

Spatial decay of G a . There is a constant k < 00 s.t. 

\\(x)^G a (x) n2 \\ < k, m+n 2 = 0,...,5, 



(2.22) 



E(m)=E 



e > s.t., for < e < e , 




p Jd (E)T e (u;,E)pj d (E)> lE p Jd (E), 



(2.23) 
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for some strictly positive constant 7^. We set 

7 := min{7£ | E G cr p (H p ) s.t. E = E(m) for some m G Jd} > 0. 

(2.24) 

Remarks. 1) 7 is exponentially small in (3, as observed in Remark 3) 
after (2.21). 

2) The operator T e (u, E) is a decreasing function of r, and an increasing 
function of R. Thus, we may assume without loss of generality that 7 is 
independent of r < 1, R > 2. 



2.1.2 Reference state a/ cf 

The reference state of the system is given by the product state 

C,/ ef = ujP ® ^ (2.25) 

where u; p is a state on B(H P ), determined by a strictly positive density matrix 
p p > 0, i.e. 

u p (A) = tr{p p A), (2.26) 

for any A G B(H P ). The state cj^ is the /5-KMS state of 21/ w.r.t. the 
free field dynamics (2.9) determined by the expectation functional (2.7). It 
describes black body radiation of the field at temperature 1/(3. 

Let (H, 7T/3, O ref ) be the GNS representation of (21, w rcf ), i.e. H is a Hilbert 
space, up is a *-morphism 21 — > B(H), and fi rcf is a vector in H s.t. ^(21)^^ 
is dense in H, and 

W rcf (A) = (tt Ici , 7T^(A)Q ref > , A G 21. 

An explicit realization of the GNS representation is well known. It was 
first constructed by Araki and Woods, [AW], and has been used recently by 
several authors. Here we just recall the explicit formulas that are useful in 
the present paper and refer to [JP], [FM] for a more detailed discussion. 
The representation Hilbert space is 

n = n p ®n p ®j=, (2.27) 
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where T is a shorthand for the Fock space 

T = T ((L 2 (R x S 2 , du x tffi)) , (2.28) 

du being the Lebesgue measure on R, and d£ the uniform measure on S* 2 . 
Here .F(^f) denotes the Bosonic Fock space over a (normed vector) space X, 

F{X) := C © (<SX 0n ) , (2.29) 

n>l 

where <S is the projection onto the symmetric subspace of the tensor product. 
We use standard notation, e.g. Q is the vacuum vector, [ip] n is the n-particle 
component of ip G J-(X), dT(A) is the second quantization of the operator 
A on X, N = dr(l) is the number operator. 

The representation map irp : 21 — > £>(7i) is the product 

where the *homomorphism 7r p : 2t p — > B(H P <S> 'Hp) is given by 

n p (A) = A® l p . (2.30) 
The representation map : 21/ — > £>(JF) is determined by 

7if (a(/i)) = / dt h(t) Tr^af (a)), (2.31) 

where 7r^j : 2U — > i3(JF) is a representation of the Weyl algebra given by 

71 w = 7T Fock O lp. 

Here, Tp is the Bogoliubov transformation, mapping 2U(Lq) to 2U(L 2 (Mx S 12 )) 
defined by W(f) i-> Wfo/), with : L 2 (M+ x S 2 ) -> L 2 (M x S 2 ) given by 



(W)(«,S) = A C^(^%f' ^ M> °' (2.32) 

v p 7V 7 Y 1 - e^" [ -V-w S), m < 0. v ' 

Remarks. 1) It is easily verified that Im (rpf, ^ / 3fi') L 2 (Rx5 2 ) = Im (/, g) L 2 (R+xS2) , 
for all /, g e L 2 ,, so the CCR (2.5) are preserved under the map rg. 
2) In the limit (3 — > 00, the r.h.s. of (2.32) tends to 

«/(«,£), M >o, , 2 33 x 

0, u < 0. I • J 
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Notice that L 2 (R + x S 2 ) ©L 2 (R + x S 2 ) is isometrically isomorphic to L 2 (R x 
S* 2 ) via the map 

(/, g) i y h, h(u, S) = ( M ^ S ^ M > °' (2.34) 

so (2.33) can be identified, via (2.34), with / e L 2 ,. Thus, 7^ reduces to the 
identity (an imbedding), tt^ becomes the Fock representation of 2H(Lq), as 
(3 — > oo, and we recover the zero temperature situation. 

It is useful to introduce the following notation. We define unitary opera- 
tors W(f) on the Hilbert space (2.27) by 

W(f) = e^ (/) , f eL 2 {Rx S 2 ), 

where <p(f) is the selfadjoint operator on T given by 

M _ m^Mi, ( , 35 ) 

and a*(f), a(f) are the creation- and annihilation operators on T , smeared 
out with /. One easily verifies that 

*Uw(f)) = W(r p f). 

The cyclic GNS vector is given by 

tt rci = tt p (g)tt, 

where Q is the vacuum in T , and 

= ^2 K<fn <S> C,plf n G H p <g> H p . (2.36) 

n>0 

Here, {/c 2 }^L is the spectrum of p p , {tp n } is an orthogonal basis of eigen- 
vectors of p p , and C p is an antilinear involution on H p . The operator C p 
comes from the identification of I 2 {Hp) (Hilbert-Schmitt operators on H p ) 
with Hp <E> Hp, via \(p){ip\ i— > </? <g> C p ^. We fix a convenient choice for C p . It 
is the antilinear involution on H p corresponding to complex conjugation of 
components of vectors in the basis in which the Hamiltonian H p is diagonal 
(i.e. it is the the time reversal operator). Then C P H P C P = H p . 
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2.1.3 Vr*-dynamical system (%Jlp,a tt x) 

Let 971/3 be the von Neumann algebra obtained by taking the weak closure of 
7173(21) in B(H), 

Wt p = B{H P ) ®\® rf(%f)" C B(H). (2.37) 

Since the density matrix p p is strictly positive, Q p is cyclic and separating 
for the von Neumann algebra n p (Ql p )" = B(H P ) <g) l p . Similarly is cyclic 
and separating for 7r^(2l/)", since it is the GNS vector of a KMS state (see 
e.g. [BRII]). Consequently Q rci is cyclic and separating for 971/3. Let J be 
the modular conjugation operator associated to (SEJt/j, f2 rcf ). It is given by 

J = J p ® J f , (2.38) 

where, for if, tjj G 7Y P , J p (</? ® = ip ® C p if, and, for *0 = {[^] n }n>o G J 7 , 

[J/^]n(«l, ■■•,««) = [^]n(-«l, • • • , -«n), for 71 > 1, 

and [J/V'Jo — [■^/V'lo G C. Clearly, Jf2 ref = f2 rcf , and one verifies that 

J p 7i p (A)J p = i p ®C p AC p , (2.39) 
Jf«itiW(f))J f = W(-e'^\(f)) = W(e'^%(f)y. (2.40) 

It is not difficult to see ([FM]) that 

a t , {n p (A)) := n p (a t>0 (A)) = e itLo np(A)e- itLo , (2.41) 

for all A G 21, where L is the selfadjoint operator on H, given by 

Lo = i/ p ®]l p -]lp®// p + dr(M), (2.42) 

commonly called the (non-interacting, standard) Liouvillian. One easily sees 
that atfl is unitarily implemented in the representation 7173 as 

M^(A)) = e< c) MA)e-^ =: a[%{^{A)), 
where the regularized Liouvillian iJf* is given by 

L[ £) = L + Xn p (V^) - XJn p (V^)J. 
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An application of the Glimm-Jaffe-Nelson Theorem (Theorem 3.1) shows 
that is essentially selfadjoint on 

p = c°° <g> ® (jr(c °°(M x s 2 )) n r ) c n, (2.43) 

where JF is the finite-particle subspace (see [FM]). Moreover, from the theo- 
rem on invariance of domains, Theorem A.l, and the Duhamel formula, one 
easily sees that 

hW* L A e) =e itL x (2.44) 

in the strong sense on H, where the Liouvillian L\ is given by 

L\ = L + XI, (2.45) 
1 = ^G a ^®\®v{Tp{g a ))-% p ®C p G a ^C p ® V {e- pu/2 Tp{g a )). 

a 

The operator L\ is essentially selfadjoint on the domain V defined in (2.43) 
and defines a *automorphism group on 97tg given by 

a t)X {A) = e ltL ^Ae~ itL \ (2.46) 

An important property of L\ is that 

e itLx J = Je itL \ for all A G R. (2.47) 

2.1.4 Characterization of the (j t A -invariant normal states on njtg 

A state w on a von Neumann algebra C B(H) is called normal iff it is 
given by a density matrix p G B(H), i.e. oo(A) = tr pA, A G 971. If r t is a 
group of homomorphisms of 971 the state is called r r invariant iff uo o r t = lo 
for all tel. In order to characterize the cr^A-mvariant normal states on 9Jlp 
it is useful to introduce the natural cone V associated to (971,9, fi ref ), which 
is defined by 

V = {AJAn rci I A G Tip} C H, (2.48) 

where the bar denotes the closure in the norm of H. The following properties 
of the natural cone are the contents of the Araki-Connes-Haagerup theorem, 
a deep result in the theory of von Neumann algebras (see e.g. [BRI]). 

Given any normal state u> on 971/3, there is a unique vector £ G V s.t. 
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oj{A) = (6 A£), for all A G 971/3. Moreover, if cui and co>2 are normal states on 
9JT/3 with corresponding vectors 6> 6 in P then 

||6-6II 2 < ||wi-w 2 || < 116-611 116 + 611- (2.49) 

The norm of a state u of 9JT/3 is given by = sup AGOT(3 

It is not difficult to see that (2.47) implies that e aL *V = V, for all A G R 
and tel. From the uniqueness of the vector representative in the natural 
cone it follows that the cr^-invariant normal states are in one-to-one corre- 
spondence with the unit vectors in the set V fl ker La, which, for A = 0, is 
given by 

V n ker Lq = V fl span|(/? m ® </? n ® f2 | m, n G M.,E{m) = (2.50) 

We will show in Theorem 2.3 that, for A 7^ 0, the cr^-invariant normal states 
are given by the subset of (2.50) determined by the modes m,n G Ai\Jd that 
do not interact with the field. 



2.1.5 A quick-reference list 

For the convenience of the reader and for future reference, we collect the def- 
initions of some important operators in a list. Generally, if n is a projection 
then we set 7f = 11 — ir. 



Pd 


projection onto the discrete subspace of H p 


Pc 


projection onto the continuous subspace of H } 


Pm 


one-dimensional projection onto the mode m < 


PJ d 


= ^meJaPm 


PJc 


spectral projection of H p onto the interval J c 


p 


= Pj d + Pj c 


p 


— p <S> p <S> 1/ is a projection on H p <8> Ti, P <E> T 


pi 


— p<S>p<S>lf 


pr 


= p®p®lf 


pO 


= p®p®tj 




projection onto kerL p 


n 


= P <S> Pa is the projection onto ker L 
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2.2 Main results 



Our main results concern the dynamical system (Wlp, c tt \), where we have de- 
fined the von Neumann algebra DJlp in (2.37), and where a t ,\ is the *automor- 
phism group (2.46) of 971/3 generated by the Liouvillian (2.45). 

The following theorem describes some properties of eigenvectors of L\ 
which, as we have seen in Subsection 2.1.4, play an important role in the 
characterization of invariant normal states. 

Theorem 2.1 (Bounds on eigenvectors). 

1) Assume that either Condition A or Condition B of Section 2.1.1 holds. 
Let N = dr(H) denote the number operator on Fock space JF(L 2 (M x 
S 2 )). Any eigenvector ip of L\ satisfies ip G V(N 1 / 2 ), for any A G R, 
and there is a constant k < oo s.t. 

\\N 1/2 i>\\ < k\X\ (2.51) 

The constant k satisfies k < k'(l + 1/(3), where k' depends on the 
interaction, but not on (3. 

2) Assume Condition A. Given any < (3 < oo, there are constants 
Xo((3) > 0, k{(3) > 0, s.t. if < |A| < \o((3), and if tp is an eigenvector 
of L\, then 

\\Po®P^\\ > Hf3)Ul (2-52) 

where Pq is the projection onto the zero eigenspace of L p , P = 1 — P , 
and Pq is the projection onto the vacuum sector in T . We have Ao(/3) > 
k'j (see (2.21) and remark 2) thereafter) and k((3) > Ary 2 , for some k 
independent of (3 and A (i.e., both constants decay exponentially in (3, 
for large (3 ). 

The proof of Theorem 2.1 is given in Section 4. Here we show that the 
bounds (2.51) and (2.52) imply that bifurcations of stationary states for the 
interacting dynamics generated by L\ from any stationary states for A = 
cannot occur. 

Let II = Pq®Pq denote the projection onto the zero eigenspace of Lq and 
set II := 11 — 11. Assume that ip is an eigenvector of L\, for some < |A| < Ao- 
Using the decomposition H" = P ® Pn + Pq. and (2.51), (2.52), we have that 

||nV|| > ||p ® Pt*l>\\ - \\p^\\ > (HP) - k\x\)U\l 
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Let ipo be an arbitrary element of ker Lq. Then 

U - ^|| > ||n0|| — H-00 — ity|| > ||n0|| - ||0o - 011, 

so 

1100- 0|| >l(k(P)-k\X\)U\\. (2.53) 

This shows that for < |A| < min(Ao, the distance between any eigen- 

vector of L\ and any eigenvector of L is greater than k(/3)/A. Combining 
(2.53) with (2.49) yields the following result. 

Theorem 2.2 (No bifurcation). Assume that the Condition A of Section 
2.1.1 holds, and that < |A| < min(Ao, t;^— ), where \q, k{(5), k are the 
constants in Theorem 2.1, 2). For any normal o- t -invariant state u° on 9Rp 
and any normal a tt \-invariant state uj x on Wlp, 

\\uj°-u x \\ > k{(3) 2 /16. (2.54) 

Our next result shows that the modes of the particle which are coupled 
to the field do not give rise to invariant states; (compare with Section 2.1.4). 

Theorem 2.3 (Instability of normal invariant states). Assume that 
Condition B of Section 2.1.1 holds. Given any < (3 < oo and any r > 0, 
there is a Ao(/3, r) > s.t. for < |A| < Ao(/3,r) ; the a \ ^-invariant normal 
states on DJlp are in one-to-one correspondence with the unit vectors in the 
set 

Vn spanjy? m <g> ip n <g> Q \ m,n G M\J d , E(m) — -E'(n)}. (2.55) 

We have Ao(/3, r) > k^ 2 r, for some k which is independent of f3,r and where 
7 is given in (2.24)- 

3 Virial theorems and the positive commuta- 
tor method 

Our proofs of Theorems 2.1 and 2.3 are based on the positive commutator 
method, which we explain in Section 3.2. In the next section we describe an 
essential ingredient of this method, the virial theorem. 
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3.1 Two abstract virial theorems 

Let H be a Hilbert space, D C Ji a core for a selfadjoint operator Y > 1, 
and X a symmetric operator on T>. We say the triple (X, Y, T>) satisfies the 
GJN (Glimm-Jaffe-Nelson) Condition, or that (X, Y, V) is a GJN-triple, if 
there is a constant k < oo, s.t. for all E V: 

< k\\Y^\\ (3.1) 

±i{(A"^,y^) - (y^,A^)} < fc(^,y^). (3.2) 

Notice that if (Xi,Y,£>) and (X 2 ,Y,£>) are GJN triples, then so is (Xi + 
X 2 ,y, V). Since y > 1, inequality (3.1) is equivalent to 

|W||<MY^||+W||, 

for some fci, k 2 < oo. For a more detailed exposition of the following results 
(and proofs of Theorems 3.2 and 3.3) we refer to [FM]. 

Theorem 3.1 (GJN commutator theorem) If (X,Y,V) satisfies the 
GJN Condition, then X determines a selfadjoint operator (again denoted by 
X), s.t. T>(X) D T>(Y). Moreover, X is essentially selfadjoint on any core 
for Y, and (3.1) is valid for all ip G V(Y). 

Based on the GJN commutator theorem, we next describe the setting for 
a general virial theorem. Suppose one is given a selfadjoint operator A > 11 
with core V C TL, and operators L, A, N, D, C n , n — 0, 1, 2, 3, all symmetric 
on V, and satisfying 

{<p,Dil>) = i{{L V ,N^)~{N^,Li,)} (3.3) 
C = L 

(<p,C n ^) = i{{C n ^<p,A^)-{A<p,C n ^)}, n= 1,2,3, (3.4) 
where </?, ip E V. We assume that 

• (X,A,X>) satisfies the GJN Condition, for X = L,N,D,C n . Conse- 
quently, all these operators determine selfadjoint operators, which we 
denote by the same letters. 

• A is selfadjoint, T> c T>(A), and e leaves T>(A) invariant. 
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Remarks. 1) From the invariance condition e itA V(A) C P(A), it follows that 
for some < k, k! < oo, and all ^ G 22(A), 

\\Ae itA ijj\\ < ke k 'W\\A^\\. (3.5) 

A proof of this can be found in [ABG], Propositions 3.2.2 and 3.2.5. 

2) Condition (3.1) is phrased equivalently as "X < kY, in the sense of Kato 
on £>" . 

3) One can show that if (A,A,TJ) satisfies conditions (3.1), (3.2), then the 
above assumption on A holds; see Theorem A.l. 

Theorem 3.2 (1 st virial theorem) Assume N and e ttA commute, for all 
t eR, in the strong sense on V, and that 

D < kN 1 / 2 , (3.6) 
Ci < kN p , for some 0<p<oo, (3.7) 
C 3 < kN 1 ' 2 (3.8) 

in the sense of Kato on V, for some k < oo. Then, if ip G T>(L) is an 
eigenvector of L, there is a family of approximating eigenvectors {ip a } C 
T>(L) n V{Ci), a > 0, such that ip a — > ip (inH) as a — > 0, and 

\im(^ a ,C^ a ) = 0. (3.9) 

Remarks. 1) It is not necessary that N and e' ltA commute for the result 
to hold, but this will be the case in all our applications (see [FM] for the case 
where N and A do not commute). 

2) In a heuristic way, we understand C\ as the commutator i[L, A] = %{LA — 
AL), and (3.9) as (tp, i[L, A]ip) = 0, which is the usual statement of the virial 
theorem; see e.g. [ABG], and [GG] for a comparison (and correction) of virial 
theorems encountered in the literature. 

The Virial Theorem is still valid if we add to the operator A a bounded 
perturbation Aq leaving the domain of L invariant. 

Theorem 3.3 (2 nd virial theorem) Suppose that we are in the situation 
of Theorem 3.2, and that A is a bounded operator on H, s.t. RanA C 
T>(L) fl RanP(A < no), for some n < oo. The commutator i[L,A ] = 
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i(LAo — AqL) is well defined in the strong sense on T>(L). For the same 
family of approximating eigenvectors as in the previous theorem, we have 
that 

lim (V> Q , (Cx + i[L, A ])4> a ) = 0. (3.10) 

3.2 Outline of the proofs of Theorems 2.1 and 2.3; the 
positive commutator method 

The positive commutator method gives a conceptually easy proof of the ab- 
sence of point spectrum of L. We outline a version that is adapted to the 
proofs of Theorems 2.1 and 2.3. The full proofs are given in Sections 4 and 
5. In the present section we use the notation of Section 3.1 and write i[L\, A] 
for d. 

Outline of the proof of Theorem 2.3. According to the discussion of 
Section 2.1.4 we have to show that PflkerL^ is given by the set (2.55). The 
Liouvillian L x is reduced by the decomposition 

H = Ran P° © Ran P © Ran P l © Ran P r , 

where the various projections are defined in Section 2.1.5. It is easy to see 
that V D RanP' = V D RanP r = {0} and that L x \ RanP = L \ RanP . 
Consequently, V D ker L\ = V fl (ker L \ Ran P° U ker L\ \ Ran P) , and to 
prove the theorem, it is enough to show that 

kerL A \ RanP = {0}. (3.11) 

We construct selfadjoint operators A, A, A such that the conditions of Sec- 
tion 3.1 are fulfilled, with L — L\ and iV = dr(ll). The operators A and Aq 
have the properties that 

i[L x , A] + i[L, A ] > PM P + Mi, 

where the bounded operators M and M 1 satisfy 

PMxP = 0, 
(ij,M ij) > 5U\\ 2 , (3.12) 
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for any ip G kerZ-A \ RanP, and where 5 is strictly positive. From Theorem 
3.3 we obtain 

= lim (V> Q , (i[L x , A] + i[L x , A \) > 5|H| 2 - 

a 

Because 5 > 0, we have that kerZ-A \ RanP = {0}. 

We now explain how to arrive at the key inequality (3.12). The form of 
the non-interacting Liouvillian, Lq, given in (2.42) suggests to consider 

A = x A P x ® \ ® If - l p ® xA P x ® 1/ + 1 P ® Hp ® dT(id u ), 

where A p is the dilation generator, see (3.40), and dT(id u ) is the second 
quantization of the translation generator in the radial variable u of L 2 (M x 
S 2 ,du x dE), see (2.28). Here, x is a function of H p with support in an 
interval in R + , containing J c but not {0}, and such that x\j c — 1- Then we 
have 

i[L x , A] = x 2 H p ® l p ® l f + l p ® x 2 H p ®l f + N + U + XI 1 , 

where N is the number operator, U is a non-negative operator, because of 
the choice of the parameter \i in the potential v of equation (2.2), and where 
l\ = i[I,A] is infinitesimally small w.r.t. N, 

±\I 1 <cN+—k, (3.13) 

c 

for any c > and some k < oo. The role of x is to project out the discrete 
modes in Jd- Using that P = p®p®lf, px 2 = pj c , and that p.j r H p > rpj c , 
since the interval J c is away from the origin by a distance of at least r, we 
obtain 

Pt[L x ,A]P 

> P (pj c H p ® lp ® P n + l p ® p Jc H p ®P n + ^Pn^j P - k\ 2 P 

> min(r, 1/2)P (p Jc ® l p ® P n + l p ® p Jc ®P n + ^P^j P - k\ 2 P, 

where we choose c = 1/2 in (3.13). The first term on the r.h.s. is strictly 
positive except on the subspace Ranp Jd ®pj d ® Pq C RanP. We decompose 

p.j d ®p Jd ®Pn = Pn+ J2 Pm®Pn®Pn, (3.14) 

m,n£ J a 
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where II is the projection onto the kernel of Lq. Note that PIT is finite- 
dimensional. On the range of the second projection on the r.h.s. of (3.14), 
the free Liouvillian satisfies 

\L \ > min{|P(m) - E{n)\ \ m,n e J d ,E(m) ^ E(n)} > 0. (3.15) 

Let A be an interval around zero whose size |A| is smaller than the r.h.s. of 
(3.14), and let E A be the spectral projection of L onto A. Then we have 
that E A pj d <S> pj d <S> Pq — E A PU = PH. Consider the decomposition 

Ran PE° A = Ran PE° A P © Ran PIT. (3.16) 

From the above discussion it is apparent that on the block Ran PE A P, 
i[L\,A] is bigger than r — k\ 2 > r/2 (we require |A| < k\/r), while the 
commutator is zero on the block Ran PIT. This is where we introduce the 
operator Aq. 

One can choose A s.t. PUi[L\, A ]UP is strictly positive provided the in- 
teraction satisfies the Fermi Golden Rule Condition. Moreover, on Ran PP^IT, 
i[L\, Aq] is small relative to r. The construction of Aq has been given, in the 
context of zero temperature systems, in [BFSS], and has been modified for 
positive temperature systems in [M]. 

The above discussion shows that the operator A] + i[L A ,A ] has 
strictly positive diagonal blocks in the decomposition (3.16). An application 
of the Feshbach method then shows that 

E° A P (i[L x , A] + i[L x , A \) PE° A =: E° A PM PE A (3.17) 

is strictly positive. Since H-E'aV' — V'll ^ ^l-^l H^ll; f° r an y tp ^ kerL^, we can 
pass from (3.17) to estimate (3.12). 

In this proof, the coupling constant cannot be chosen idependently of the 
inverse temperature (3. This is due to the fact that the constant 5 in (3.12) 
is proportional to 7, see (2.24), which in turn decays exponentially in (3. We 
have to require that certain error terms which depend on A are small w.r.t. 
7, hence the /5-dependent smallness condition on A. 

Outline of the proof of Theorem 2.1 . Part 1) is an easy consequence of 
the virial theorem combined with the bound 

i[L x ,A] > l -N-k\\ 
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for A = o\T(id u ). The proof of part 2) proceeds as follows. We construct Aq 
(the same as for the proof of Theorem 2.1) s.t. 

i[L x , A] + i[L x , A ] > Kl U - k 2 P <g> P n , (3.18) 

for some K\,k 2 > 0, and where II is the projection onto the kernel of 
L . If ip is an eigenvector of L\ then by part 1) we have that ||IT/>|| > 
(1 — fclADUV'll — ||Po ® PfiV'll- Inserting this bound into (3.18) and using the 
virial theorem yields the bound (2.52). 

This outline indicates that the proofs of Theorems 2.1 and 2.3 consist of 
two steps. First we verify that the virial theorems are applicable and second 
establish a positive commutator estimate in the above sense. The latter task 
is carried out in Sections 4 and 5. 

3.3 Applications of the virial theorems 

Corresponding to the different hypotheses of Theorems 2.1 and 2.3, we intro- 
duce two sets of operators A, L, A, A , N, and verify, in each case, that the 
virial theorems are applicable. The following objects appear in both applica- 
tions: the Hilbert space is the GNS space given in (2.27); the dense domain 
V is chosen to be 

V = C^°(R 3 ) <g> C °°(K 3 ) <g> V f , (3.19) 

where 

V f =F(C?(R x S 2 )) njFo, 

where the Fock space T has been defined in (2.29), and JF denotes the finite- 
particle subspace. The operator L is the interacting Liouvillian introduced in 
(2.45), and N = dr(ll) is the particle number operator in T = JF(L 2 (Rx S 2 )). 
Clearly, X = L,N are symmetric operators on T>. The operator D, defined 
in (3.3), is given by 

a 

-\ ® C p G a , # C p ® (- a *(e-M 2 T P (g a )) + a( e -^\(g a ))) }. (3.20) 

We define a bounded, selfadjoint operator A on 7i by 

A = ie\(IlIR 2 Ti-TlR 2 m), (3.21) 
R\ = (Ll + e 2 )- 1 . (3.22) 
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Here, 9 and e are positive parameters, and II is the projection 



IT = P ®P n , 
= P(L p = 0), 

TT = i-n. 



(3.23) 
(3.24) 
(3.25) 



We also introduce the notation 



R e = UR e . 



Notice that the operator Aq satisfies the conditions given in Theorem 3.3 
with n — 1. Moreover, [LA ] = LA — A L extends to a bounded operator 
on the entire Hilbert space, and 



This choice for the operator A was initially introduced in [BFSS] for the 
spectral analysis of Pauli-Fierz Hamiltonians (zero temperature systems), 
and was adopted in [M] to show return to equilibrium (positive temperature 
systems). The key feature of A is that 



is a non-negative operator. The Fermi Golden Rule Condition, (2.20) (or 
(2.23)), says that it is a strictly positive operator on Ran II. 

Proposition 3.2 Assume (2.20) and let < e < e . Then 



where 7 is given by (2.21). Assuming condition (2.23) instead of (2.20), 
the same lower bound holds (with 7 given in (2.24)) if we replace II by IIP 
(P — p <S> p, p — pj d + pj c ) and I by the regularized interaction; see (2.13). 

The proof of Proposition 3.2 is given in Section A. 4. 
Next, we define the operators A and A and verify the hypotheses used in 
Section 3.1. 




(3.26) 



iU[L,A ]U = 29\ 2 UIR 2 JU 



uiR 2 e m > -711, 



(3.27) 
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3.3.1 Setting for Theorem 2.1 

We define 

A = A p ®l p ® l f + l p ® A p ® l f + l p ® J,,® A f , (3.28) 

A p = H P -E Q + 1, (3.29) 

A f = dr(w 2 + l), (3.30) 

where, we recall, E = inf <j(H p ) < 0. Clearly, A is essentially selfadjoint on 
the domain V defined in (3.19), and A p > 1, Af > 1. In what follows we 
shall often use the standard fact that if / e L 2 (R x S* 2 , du x dE), then a # (/) 
is relatively N 1 ^ 2 bounded in the sense of Kato. This implies immediately 
that a*(/) is relatively A 1 ^ 2 bounded. 

We verify that (L, A, T>) is a GJN triple. The bound (3.1) is trivial by 
the above observation, and the fact that Tp(g a ) € L 2 (R x S 2 ). Next, the only 
contribution to the commutator of L with A comes from the interaction, and 
a typical term to estimate is of the form [G a , H p \ (g l<g ip(Tp(g a )) + G a <g> H p <g> 
[</?(Tg((? a )), A/]. Using the bound (2.18), we obtain for the first term 

|(V>, [C a , ifp] (g) l p (g) ip{Tp{g a ))ij})\ 

< k\\[G a ,H p ](H p -E + l)-V 2 \\ \\A l J 2 ®^|| ||AyV|| 

< k{ip,Aip). (3.31) 

Next, 

|^,G Q ®]1 P ® [^(^(^)),A / ]^)| < fclK^ + i) 1 / 2 ^)!!^^) IIA^VH 2 

< k{ip,Aii>), (3.32) 

where we have used that 

[a'fofoO)^] = a* ((« 2 + 1)^)) , 
[a(r^(c/ a )), A/] = -a ((u 2 + 1)^^)) , 

so that [<p(Tp(g a )), Af] is still N 1 ^ 2 bounded, since Tp(g a ) has the decay prop- 
erty (2.17). The form bound (3.2) follows from these observations. In a 
similar way, one shows that (D,A, V) is a GJN triple. 
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Next, we define the operator A = Aj to be the selfadjoint generator of 
the translation group acting on the radial variable of elements in T by 

[e uAf i)} n {ui, Ei, . . . , Un, S n ) = [ip] n (ui - t, Ei, . . . , - t, E n ), i G R. 

In what follows, we will often not display the angular variables Ei, . . . , E n . 
We set e ltA fQ := f2. Clearly, "D C V{Af), V is invariant under e ltA ! , hence a 
core for Aj, and Aj acts on D as 

A/ = dF(id u ). (3.33) 

An easy calculation shows that, on V, 

Ae itA f = e tA s ^ A + dT ( 2ut _ t 2^ ^ ^ 3 3 ^ 

so estimate (3.5), with k' = 0, is satisfied for all G X>, hence for all ^ G 
P(A). For A := A/, we find that 

d = N + XI U (3.35) 

C 2 = A/ 2 , (3.36) 

C 3 = AJ 3 , (3.37) 

where 

I n = i n J2 i Ga J P ¥>(H0„) n T/j(0a)) 

-Hp ® C p G a C p ® v^((-^) n e^ u/2 r /3 (^))}. (3.38) 

We now show that (C n , A, X>) are GJN triples, for n = 1, 2, 3. The opera- 
tors I n are A 1 / 2 -bounded, since Tp{g a ) and e _/3u / 2 Tg((? Q ) are in in the domain 
of the operators (id u ) n , n = 1,2,3 (see also (2.16)), hence (3.1) holds. Note 
that this also yields (3.7). Next, we need to calculate the commutators of C n 
with A. The estimates on the commutators of I n with A, for n — 2, 3, are 
similar to the ones for n — 1. The latter has been outlined above; it requires 
that (u 2 + 1)(— id u ) n T[3(g a ) G L 2 (R x S 2 ), which is guaranteed by conditions 
(2.16) and (2.17). 

This discussion shows that we are in the situation described in Section 
3.1, and Theorems 3.2, 3.3 apply. 



29 



3.3.2 Setting for Theorem 2.3 

We define the operator A as in (3.28), but where A p is now given by 

A p = -A + x 2 . (3.39) 

A is essentially selfadjoint on V (see (3.19)), A p > 1, A/ > 11. 

Verifying that (Lj,A,V) is a GJN triple is done as in Subsection 3.3.1, 
using that [G a ,.j, A p ] is bounded; see Lemma A. 5. It is also easy to check 
that (D, A, V) is a GJN triple. 

Next, we define an operator A differing substantially from the choice 
A = Af (see (3.33)) in Subsection 3.3.1: We add a (regularized) dilatation 
on the particle space to Af. 

Let x be a smooth characteristic function of the set J c (with 

the property x\j c — 1), which has compact support not containing zero. We 
define x{H p ) = f x( s ) e%sHp > where x is the Fourier transform of Xi an d we 
abbreviate x{H P ) by x- Let A p be the symmetric operator on C£°(IR 3 ) given 
by 

A p = - l -{x- V + V-x). (3.40) 

Notice that (A p , A p , C^(R 3 )) is a GJN triple, so A p is essentially selfadjoint 
on C^°(IR 3 ). We denote the selfadjoint closure again by A p . 

Remark. To show that A p is essentially selfadjoint on C^°(IR 3 ), we can 
also use the fact that the dense set C^°(1R 3 ) is invariant under the group of 
dilatations on L 2 (IR 3 , d 3 x), hence a core for the selfadjoint generator of this 
group. The generator acts on C£° as in (3.40). 

Proposition 3.3 (xA P x, A p , C^°(IR 3 )) is a GJN triple. In particular, xA p x 
is well defined and symmetric on C^°(R 3 ), and it is essentially selfadjoint on 
C£°(R 3 ). We denote the selfadjoint closure again by xA P X- 

We give the proof in Section A. 2. Let us now define the operator 

A = X A P x ® H P - Hp ® XAX + A f , (3.41) 

which is essentially selfadjoint on T>. It follows immediately from Proposition 
3.3, Theorem A.l, and relation (3.34), that e ltA leaves V(A) invariant, and 
that the estimate (3.5) holds true. We calculate explicitly 

C n = S nA N + adg pX (# p ) <g) l p + (-l) n l P ® adg pX (# p ) + XI n , (3.42) 
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for n = 1,2,3, where we define the multiple commutators ad y (X) = X, 
and for n > 1, ady\x) = i[ady _1 \x), Y], in the weak sense on C£°(1R 3 ) x 
C^°(]R 3 ). For n = 1, 2, 3, we have defined 

^ = E (f) 2 " (n " fe) E { ad xXx(^) ® ^ ® v 

+ (-1)-% ® adj£?(G«) ® V? (H^e^V^)) } . (3.43) 

Note that 

ad^ x (flp)=x(^ + W) X , (3.44) 

with 

^ = 4 (I .V. + 2.) = I(^ + (1-^), (3.45) 

and the choice of (/,, p given in (2.2) implies that 

W > 0. (3.46) 
In Appendix A. 3 we prove the following Proposition. 

Proposition 3.4 (C n ,A, V) are GJN triples, for n = 1,2,3, and the esti- 
mates (3.6) -(3.8) are satisfied. 

This shows that with the choice of operators introduced in this section, The- 
orems 3.2 and 3.3 apply. 



4 Proof of Theorem 2.1 

1) Set h = i[I,A f ], where A f and / are given in (3.33), (2.45). For -ip e 
P(A 1 / 2 ), we have 



< 



+ 



< 2\\I 1 N- 1/2 P n 



\N 1/2 tP\\. 



This shows that in the sense of quadratic forms on T>(N), A/i > — cN — — , 
for any c > 0, where k = WhN'^Pnf < k' Z a \\d u rp(g a )\\ 2 L2 < k'(l + l/(3). 
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Let C\ = i[L\,Af] = N + Ap and choose c = 1/2. Then we find C\ > 
|7V — /cA 2 , in the sense of forms on T>(N) C P(Ci), and from Theorem 3.2 

= lim U a ,C^ a ) > ^ Hm ||iV^J 2 - fcA^f, 
where ip is an eigenvector of L\, and ?/> Q its regularization. It follows that 

lim || N l/ Va II 2 < A;A 2 ||^|| 2 , 

which tells us that ^ e V^N 1 / 2 ), and that ||iV 1/2 V;|| < fc|A| ||^||. 

2) In what follows, the constants k, hi, Ai, A2 are independent of /3 > /?o, 
where /5o > is arbitrary but fixed. In the course of the proof we will impose 
several conditions on the parameters e, A, 9 which are collected in (4.14). We 
adopt the notation of Subsection 3.3.1. 

On V(N) C V(Ci), we define the operator 

B = C 1 + i[L x ,A ]. 

Recall that A is defined in (3.21), that we write R £ = UR e , and that 
P n IP n = PnhPn = 0. Using that P n = U + P <g> P a , one finds that 

P n BP n = UBU + P ®P n BIl + ILB 'J 5 ®P n + P ®P n BP ®P n 

= 29\ 2 TlIR 2 e m + 9X 2 (P <g> P Q IR 2 e IU + UIR 2 e IP <g> P n ), (4.1) 

P^PPn = \PnhPn + 9\UIR 2 € L P n + 9\ 2 UIR 2 JP n , (4.2) 
P^PPn = XPnhPn + 9\P n L R 2 m + 9\ 2 P n IR 2 JU, (4.3) 

p n pp n = p n N + p n (a/i - 0A 2 (/n/p 2 + p 2 /n/)) p q . 

From the estimates ||^ A^" 1 / 2 /! A^" 1 / 2 ^ || < k, \\IUI\\ < k, ||P 2 || < e~ 2 , we 
see that there is some constant Ai < 00 (independent of A, e, 9), s.t. 

PnBPn > ^Pq, (4.4) 

provided 

|A|, — <Ai, (4.5) 
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see also (4.14). Using the estimates 

\\IR 2 e IUl \\UIR 2 £ I\\ < e~ 2 k and \\P n h\\, \\hPn\\ < k, 

where k is independent of the parameters A, 9, e, we arrive at the following 
lower bound. For any G V(N) and some k\ < oo 

(B)+ > 26\ 2 (mR 2 m)^+ 1 -\\p n( f ) \\ 2 

a\2 

-fci — ||IT0|| ||P ®Pn^||-fci|A| \\Pq4W 
^( 2 e\\\+2k 1 —) \\R e m<f>\\ ||P n 0||. (4.6) 



Clearly, ||P e /n0|| \\P n <j>\\ < 5 (UIRJUJ + S^WP^W 2 , for any 5 > 0. 

Choosing appropriate values of 5, we bound the last line in (4.6) from below 
by 

2 / TT T ID 2 TT1 \ A (a , L2^^ MO i||2 



#A 2 (n/p;/n) ~a(o + k 2 ^ ) ||p. 



and it follows that 

(B), > ^7l|n0|| 2 +Q-4^-4A;^||P,0|P 



d\2 _ _ 

-h — WW \\Po ® Pn<f>\\ - \\PnH \\PnH, (4.7) 

where we have used (3.27) and hence assumed that < e < eo- Using that 
< I|n0|| + ||Po ® Pq0||, we estimate the two terms in the last line on 
the r.h.s. of (4.7) as 

-fcxlAI ||Pn0|| ||P^|| > -i||P^|| 2 -8A 2 A; 2 (||n0|| 2 +||P o ®P n 0|| 2 ), 

n\2 _ i n\2 n\2 

-h — IIU0H ||P o ®P n 0|| > --— 7 ||n0|| 2 -2A; 2 — 7 ||P o ®P f ,0|| 2 . 
Using these two estimates in (4.7), we arrive at 

1 & ^ t 1 1 1 1 . ,iiO ^ ■» O t o / , ^ 



(B), > A 2 ( —7 - 8A; 2 j ||n0|| 2 - 2A 2 A; 2 (4 + ^ ) ||P ® P,0|| 2 , (4.8) 
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where we require the condition 

]-49-4kf-^->0, (4.9) 
4 e z 

which guarantees that the contribution of the term in (4.7) which is propor- 
tional to ||Pn0|| 2 is non-negative, and can hence be dropped. (4.9) is satisfied 
if (4.14) holds. Let ^ = ^6 V(N) be the regularization of the eigenvector 
ip as defined in Theorem 3.2. Then it follows from (4.8) that 



= lim(£), ^ilin^f-^HPo®^!! 2 , (4.10) 



where 



«i = l -- el -%kl>kl (4.11) 

k 2 = 2k\ (^4+ -f^J > 0. (4.12) 

The lower bound (4.11) is a consequence of e < yfp-, see (4.14). 
From (2.51), we find that 

||Ity|| > |H| - \\P n i/>\\ - ||P <E> Pn^\\ > (1 - fc|A|)||^|| - ||P ® Pn4>\\. 

Thus there is a positive constant A2 (independent of e, A, 9) s.t. if < |A| < A2 
then ||Ity|| > - \\P <g> P n ip\\. Thus we get from (4.10) 

I|7oWII>^ 1 + (k ^ i)1/2 IM|. (4.13) 
Consequently, under the conditions that 

1 <- • l±L 

VeV" 1 ' lyftkjy ~~ 32' e<mm \l8A;f 



< |A| < min { Ai, A 2 , — ( V ^1 + , 7= ) } , # < ™> e < min { 7777^2 , e o 



(4.14) 

Pni>\\ > \ . 1 = (4.15) 

2 l + A /2(4 + -^) 



we obtain 



Choose for instance 9 = 1/100, e = min{ 200 7 ofc 2 ,eo}- Then (4.14) holds pro- 
vided < |A| < &7, for some k independent of (3 provided that j3 > j3o, with 
Po > arbitrary but fixed. For large f3, the r.h.s. of (4.15) behaves like 7 2 . 
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5 Proof of Theorem 2.3 



The <7 ti A-invariant normal states on DJlp are in one-to-one correspondence 
with the normalized vectors in the span of V n ker La (see Section 2.1.4). 
Our task is to show that V Piker La equals the set (2.55). 

In this section, we will always deal with the cutoff interaction (determined 
by G a j) but we shall drop the subscript j in the notation. 

5.1 Reduction of the Liouvillian 

We define the projection p = p Jd + p Jc , where pj d ,pj c are the projections 
corresponding to the discrete and continuous modes in J d and J c , respec- 
tively; (see also (2.14)). Setting p = l p — p, P = p <g> p <g> If, we decompose 
P = 1 - P as P = P l + P r + P°, where 

P l = pig) p ® tf P r =p®p®l f , P°=p®p®l f . (5.1) 

It is easy to verify that the (regularized) Liouvillian La, defined in (2.45), is 
reduced by the decomposition 

H = Ran P © Ran P l © Ran P r © Ran P°, 

and that 

L A \ RanP = L \ RanP . (5.2) 
From the definition, (2.38), of the modular conjugation J, it follows that 

jpl = prj ^ 5 _ 3 ) 

Because every ip e V satisfies Jip = ip, (5.3) implies that V fl RanP' = 
V Pl RanP r = {0}, and, consequently, we have that 

V n ker L A = V n (ker L A \ Ran P° U ker L A f Ran P) . (5.4) 

We prove in the next section that kerL A f RanP = {0}, which, together 
with (5.4) and (5.2), shows that V nkerL A is given by the subspace defined 
in (2.55). 



35 



5.2 The kernel of L\ \ RanP 

Theorem 5.1 Given any < (3 < oo and any r > there is a \ ([3,r) > 
s.t. ifO< |A| < A (/3,r) i/ien kerL A f RanP = {0}. Here, X (P,r) > k'fr, 
for some k independent of (3,r. The constant 7 is given in (2.24) ■ 

Proof. We use the notation of Subsection 3.3.2 and write L for L\. In the 
spirit of the positive commutator method outlined in Section 3.2, we want to 
establish a lower bound on the expectation value (ip,(C 1 + i[L, A ])ip) (see 
also (3.42)), where ip is a (hypothetical) eigenvector of L in RanP. 

Using the relative bound 

A 2 

±A/i > -cN -—, Vc> 0, (5.5) 

with c = 1/10, and (3.46), we obtain a lower bound (always in the sense of 
quadratic forms on V) 

Ci + i[L,Ao] 

9 A 2 

> x 2 H p ® Hp + Hp ® x 2 H p + —N + i[L, A ] - — 

( 9— A 2 A 

> P [x 2 H p ® \ ® Pn + \ <g) x 2 ^ P ® Pn + — Pn + i[L, A ]-—jP 

A 2 — 

+P^[L, A ]P + Pi[L, A ]P + Pi[L, A ]P - —P 
=: PM P + M U (5.6) 

where the bounded operators M Q and Mi are given by 

9 — 

M := p Jc H p ®l P ®P {l + l P ®p Jc H p ®P n + —P n 

+i[L,A ]-^, (5.7) 

M x := Pi[L,A }P + Pi[L,A }P + Pi[L : A }P-^P. (5.8) 

The difficult part of the proof of Theorem 2.3 is contained in the following 
two propositions. 



36 



Proposition 5.1 Suppose Proposition 3.2 holds and that the parameters sat- 
isfy 

< | A| < min ^1, y/r, ~^=j , < e < min(56>7, e ), < 9 < r/32, 

3 (5-9) 
where k is a constant depending on the interaction, but not on any of the 
parameters e, A, 6, nor on {3 (for (5 > with j3 > fixed). Then there is 
an interval A around zero such that 

OX 2 

PE° A M E° A P > —^E° A P, (5.10) 

where E A = E A (L ) is the spectral projection of L onto A. 

Proposition 5.2 Assume that the conditions of Proposition 5.1 are satisfied 
and that 

|A| < Jmin(l,e,0/e). (5.11) 
If ip G Ran P is an eigenvector of L\ then 

1 f)\ 2 

<^M ^>>-— 7 |H| 2 . (5.12) 

We may choose the parameters as A = A7 2 , e = ey, with A, ? and 9 
independent of the inverse temperature (3. Conditions (5.9) and (5.11) are 
satisfied provided < |A| < kr, for some k independent of (5 and r. 

Theorem 5.1 is now proven as follows. Assume that ip G RanP is an 
eigenvector of L\, and let ip a be the family of approximate eigenvectors given 
in Theorem 3.3. From (5.6), (5.12), and using that (ip,Miif)) = 0, we obtain 

1 OX 2 

= lim (ij a , (d + i[L, A })^j a ) > -— 7 |H| 2 , 
a z e 

which is a contradiction, because the r.h.s. is strictly positive. 
Proof of Proposition 5.1. Pick A such that 

|A| < ^min - E(n)\ \ m,n G J d ,E(m) ^ (5.13) 
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The Hilbert space Ran E A P has the decomposition 

Ran P^P = RanPII © Ran P^ PIT, (5.14) 

where, we recall, II = P ® Pn is the projection onto the kernel of L . We 
analyze the spectrum of the operator PE A M E A P on Ran P^P using the 
Feshbach method. For details on the Feshbach method, we refer to [BFS]. 
Let m be a number in the resolvent set of IIPP^MoP^Pn (viewed as an 
operator on RanP^PII). The Feshbach map Pn, m applied to the operator 
PE A M E° A P is defined as 

F n<m (PE A M E A P) (5.15) 

= n (pm q p - pm TIpe a (TTpp^m p^pTT - m) _1 E A PUM P^ it, 

and has the following property of isospectrality: Let a and p denote the 
spectrum and the resolvent set of an operator. Then 

z G a(PE° A M E° A P) n p(TlPE A M E A PTl) (5.16) 
z G a(F^ m {PE A M E° A P)) fl p(TlPE A M E A PTi). 

The point of the Feshbach method is that it can be easier to analyze the 
spectrum of the operator (5.15) than the one of PE A M E A P, because the 
operator (5.15) acts on the smaller space RanPII. 

Let us examine the diagonal blocks of PE a MqE a P in the decompositon 
(5.14). It is readily verified that 

-9 A 2 

UPM PU = 29\ 2 UPIRJPU PIT. (5.17) 

10 y J 

Because of (5.13), P^Po P.j d ®PJ d ® Pn = 0. Hence 

E° A PIT = E° A PP n + E° A P {p Jd <g> p Jc + p Jr <g> p Jd + p Jc <g> p Jc ) <g> P n . 

Set Q\ = E a PPq and let Q 2 be the other projection on the right side. 
PE A M E A P is diagonal in this decomposition of Ran E A PIT. We have the 
estimates 

/ n \2 0\ 2 \ ( A 2 \ 

QiMoQ, >(----k—)Q 1 , Q 2 M Q 2 > r - — Q 2 , 



10 10 e 2 / " ' " u " - V 10 
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where k = 2||iTI/|| and we used that pj c H p > rpj c . Consequently, we obtain 
the lower bound 

/ g \2 Q\2 \2\ _ r - 

PUE A M E° A UP > min I — - — - fc— , r - - E° A UP > -E° A UP, 

(5.18) 

due to (5.9). It follows that any m < r/A is in the resolvent set of the 
operator PUE° A M E A TiP and 

(PUE A M E A UP - m)' 1 <4/r. (5.19) 

We show now that 

F n , m (PE° A M E A P) > — 7 E° A P, (5.20) 

uniformly in m < r/A, provided (5.9) is satisfied, and where the Feshbach 
map was introduced in (5.15). The bound (5.20) and the isospectrality prop- 
erty (5.16) of the Feshbach map imply (5.10). 

We complete the proof of the proposition by showing (5.20). From (5.19) 
it follows that for any ip and m < r/A 

ip, UPMoU.PE° A (UPE A M E A PU - m) ^ E A PUM PU^ 

< -||TTP^M Pn^|| 2 

r 

< ^-(1 + k\ 2 /e 2 ) U, PIUlfjTlpA , (5.21) 



r 



where we estimate T1PE A M TIP = 6\nPE A LR 2 e ITlP as 

\\UPE A M UPiP\\ < 9\X\(1 + A;|A|/e)||P e /nP^||, 

with k = \\IP(N < 2)||. Taking into account (5.17) and Proposition 3.2, we 
obtain the estimate 

F u , m (PE A M E A P)>2^ 1 ^l_^(l + A;A 2 A 2 )JnP-^nP. (5.22) 

The bound (5.20) follows from (5.22) and conditions (5.9). This finishes the 
proof of Proposition 5.1. 
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Proof of Proposition 5.2. Let < g < 1 be a smooth function with 
support in the interval A, s.t. g = 1 on the interval (— |A|/4, |A|/4), and set 
g = g(L ). Since the interaction / is relatively iV^-bounded, it follows in 
a standard way (by using e.g. the functional calculus presented in Appendix 
A) that 

11(1-50)^11 < fc|A| \\(N + 1) 1/2 1>\\ <k\X\ (5.23) 
||(1 -go)P a 1>\\ < k\X\ \\N l ' 2 P^\\ <k\\\ 2 Ul (5.24) 

where we used (2.51) in the last line. Notice also that 11(1 — go) = and that 
on Ran(l — g ) we have \L \ > |A|/4 hence ||(1 — go)Rt\\ < 4/|A|. These 
estimates are used below without explicit mention. We decompose 

PM P = g PM Pg 

+2Re(l - g )PM Pg (5.25) 
+ (l-g )PM P(l-g ). (5.26) 

Proposition 5.1 yields the bound 

(^,g PM Pg^) > ||^f = ^ 7 (1 - A;|A|) 2 |H| 2 . (5.27) 
We estimate 

2Re(l - g )PM Pg > 2Re(l - g Q )Pi[L, A ]Pg - A 2 (l - g )g P, (5.28) 
and since 

(1 - g )i[L, A ]g = -9X(1 - g ) (xIHIR] - LlfjTl - XR 2 e mij g 
we conclude that 

2Re Oj, (1 - g )PM Pg ^) > -k— + l^) |H| 2 . (5.29) 



Next, we have that 

(1 - g )M (l - g ) > -9\ 2 (1 - g ) [mm] + R 2 e IUl) (1 - <? ) - ^(1 - 9o)\ 
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from which it follows that 

(V, (1 - golWfl - goM > -k^j- (^A 4 e + |M| 2 . (5.30) 
Collecting the bounds (5.27), (5.29) and (5.30) we obtain 

(♦,™.»>>^ (i-*|a|-^ + ^))|MI 2 , 

and (5.12) follows from the conditions (5.11). This completes the proof of 
Proposition 5.2 and of Theorem 5.1. ■ 



A Appendix 

A.l Invar iance of domains, commutator expansion 

The following two theorems are useful in our analysis. 

Theorem A.l (invariance of domain, [F]) Suppose (X,Y,V) satisfies 
the GJN Condition, (3.1), (3.2). Then the unitary group generated by the 
s elf adjoint X , e ttx , leaves V(Y) invariant, and we have the estimate 

\\Ye ltx ^\\ < e fe l*l||yV||, (A.l) 
for some k > 0, and all ip G T>(Y). 

Theorem A. 2 (commutator expansion, [F]) Suppose V is a core for 
the self adjoint Y > 1. Let X, Z be two symmetric operators on T>, and define 
the symmetric operators ad^(Z) on V by 

ad^(Z) = Z, 

(^M { x\z)i) = ^{(adf- 1) (Z)^X^)-(x^,ad5- 1) (Z)^)}, 

for all i\) e T>, n — 1, . . . ,M. We suppose that the triples (ad^(Z), Y, V), 
n = 0, 1, . . . , M, satisfy the GJN Condition (3.1), (3.2), and that X is self- 
adjoint, V C V{X), e ltx leaves T>(Y) invariant, and that(A.l) holds. The 
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we have on T>(Y) 

e itx Ze -ux = z _J2 l - a dP(Z) 

n=l n ' 

- f dh ■ ■ ■ /* M_1 dt M e ltMX \z)e- ltMX . (A.2) 
Jo Jo 

The following Lemma is a consequence of the above two theorems. 

Lemma A. 1 Suppose (X,Y,V) and (ad^ (Y ) , Y, V) are GJN triples, for 
n = 1,...,M, some M > 1. Moreover, assume that in the sense of Kato 
on V(Y): iad^f ] (Y) < kX , for some k > 0. For x e S{R), a smooth 
function of rapid decrease, define x(X) = f x(s)e lsX , where x is the Fourier 
transform of x- Then leaves V{Y) invariant. 

Proof. For R > 0, set X «P0 = fix(s)e isX , then X r(X) -> xPO 
in operator norm, as i? — > oo. From the invariance of domain theorem, 
we see that Xr(X) leaves V(Y) invariant. Let ip £ T)(Y), then using the 
commutator expansion theorem above, we have 

f R 

= XR{X)Yi) + / x{s)e lsX (e- lsX Ye lsX -Y)ij} 

J-R 

rR /M—l ,_ , n 

= X r(X)Y^ - / X (s)e- X J2 M^ ad x (X) 

J ~ R \n=l U - 

+ (-l) M jT d Sl ■ ■ ■ ds M e- iSMX sd { x\Y)e iSMxS j V (A.3) 

The integrand of the s- integration in (A.3) is bounded in norm by 

k(\ s \ M + 1) (\\y^\\ + \\x^\\) < k(\s\ M + i)||i^||, 

where we used that ||ad^f ) (Y)e iSMX ip\\ < k\\Xe iSMX ip\\ = k\\Xip\\. Since x 
is of rapid decrease, it can be integrated against any power of \s\, and we 
conclude that the r.h.s. of (A.3) has a limit as R — > oo. Since Y is a closed 
operator it follows that x{X)i> £ T>(Y). ■ 
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Lemma A. 2 Let x G C£°(R 3 ), X = F 2 > 0. Suppose (X,Y,V) satisfies 
the GJN condition and define F(X),x(X) via the Fourier transform as in 
Lemma A.l. Suppose F(X) leaves V(Y) invariant. Let Z be a symmet- 
ric operator on T>, s.t. for some M > 1, and n = 0, 1, ... , M, the triples 
(ad^ - (Z), Y, V) satisfy the GJN condition. Moreover, assume that the multi- 
commutators for n = 1, . . . , M — 1 are bounded, and the M-th multicommu- 
tator is relatively X-bounded in the sense of Kato on V: there is a k < oo 
s.t. W> G V: 

||ad?(Z)V|| < fc|M|, n=l,...,M-l, 
\\^\z)n < k\\X1>\\- 

Then the commutator [x(X), Z] = x(X)Z — Zx(X) is well defined on V and 
bounded: there is a k < oo s.t. \\[x(X), Z]ip\\ < k\\ip\\, Wip G V. 



Proof. We write F,x instead of F(X),x{X). Since F leaves T>(Y) 
invariant we have in the strong sense on T>(Y): 

[X,Z]=F[F,Z] + [F,Z]F. 

We expand the commutator 

[F, Z] = J F(s)e isX (Z - e- lsX Ze isX ) 

^ n=l 

dsi • • • J ds M e~ lSMX ^ ] {Z)e tSMX | . 

Multiplying this from the left with F (and noticing that F commutes with 
e lSMX ) ) we see immediately that [F, Z]F is bounded. Next, for any <p,if) in 
the dense set T>, we have the estimate \{<j), F[F, Z]ip)\ = \([F, Z]F<f>, ip)\ < 
kU\\ hence ||F[F, Z}^\\ = sup ^ eV UW' 1 ] (0, F[F, Z\f) \ < fc||^||. ■ 

A. 2 Proof of Proposition 3.3 

Before proving Proposition 3.3, we show certain triples satisfy the GJN con- 
ditions. 
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Lemma A. 3 The following triples satisfy the GJN Condition (3.1), (3.2): 

(H P ,A P ,C™(R 3 )) (A.4) 
(A p , A p , C^(R 3 )) (A.5) 
(adP p (Ap),A p ,C™(R 3 )), n=l,2. (A.6) 

Proof. For ijj G C^°(IR 3 ) we have \\H p if;\\ 2 < 2||(-A)^|| 2 + 2||t»||^||^|| 2 . 
The first term on the r.h.s. is bounded from above by 2||A p -?/>|| 2 +4Re (0, Ax 2 ijj) 
and we have the estimate Re (0, Ax 2 ip) = X^ =1 ((0> XjAxjip) + (-0, djXjip)) < 
Ej 3 =i ll^ll \\xji/>\\ < (ip,(-A + x 2 )^). Therefore, since A p > 1, it follows 
that that ||i? p ^|| 2 < /c||A p -0|| 2 . In a similar way it is simple to verify that 
±ad^ (Hp) < kA p . This shows that (A.4) satisfies the GJN conditions. The 
proof for (A.5) is similarly easy. 

From the above calculations, and ad^(A p ) = — a,d^(H p ), we see that 

||ad^(A p )^|| < fc||Ap/ 2 ^||. Moreover, we calculate 

ad« (ad«(A p )) 

= - (2x ■ Wv + 2d m v mn d n + d n v nmm + 4 A + 4x 2 ) + h.c, 

where Vki m = dkdid m v, etc. Since all the derivatives of v involved are 
bounded, the r.h.s. is again relatively A p -bounded, in the form sense on 
Co°(lR 3 ). This shows (A.6) for n = I. For n = 2, we calculate 

adg(A p ) = - (Vv ■ Vv + 2x ■ Vv + 2d m v nm d n + d m v mnn + 4A)+h.c. (A.7) 

Notice that ad^(A p ) is relatively i/ p -bounded, hence relatively A p -bounded, 
in the sense of Kato on C^°(IR 3 ), because Vt> • Vi> and x ■ Vt> are bounded. 
Moreover, one can calculate ad^ ^ad^(A p )j, and see that it is A p form 
bounded on C^°(1R 3 ). This shows (A.6) for n = 2. ■ 

Proof of Proposition 3.3. We start by showing that x leaves T>(A P ) 
invariant. This follows from Lemma A.l and the following two facts: firstly, 
(H p , Ap, C£°(R 3 )) and (adjg, A p , C^°(1R 3 )) satisfy the GJN Condition, for n = 

1, 2 (see Lemma A. 3), and secondly, ±ad^(A p ) < kH p , in the sense of Kato 
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on T>(A P ) (see (A. 7)). This shows that x '■ 'D(Ap) — > V(A P ), and consequently, 
xA P x is well defined on V{A P ). 

We have in the strong sense on T>(A p ) 

X A pX = X 2 A p + x[A P ,x]- (A.8) 
Let us estimate each term on the r.h.s. separately. For ip e C^°(M 3 ), 

||X 2 ^|| = || X 2 (V -x + 3/2)^|| < £ Wx 2 d n \\ ||^|| + ^\\ 

n 

< kj2(^, ^> 1/2 + lu\\<k(ii, (-a + x 2 )^> 1/2 , 

n 

where we used that ||x 2, 9ri|| < °°. Consequently, 

||x 2 ^||<A;||AyVl|. (A.9) 

Next we have on T>(A P ) 

x\A p ,x] = ~xj ds X (s) jy Sl e^ H ^d { }l(H p )e'^ H ^ 

where ad^(if p ) = H p + W, see (3.44), and since x commutes with e lSlHp , we 
obtain the estimate 

||X[A,, X ]||< ^3 + ||W||oo)<c30. (A.10) 

It follows, together with with (A.8) and (A.9), that 

\\ X A pX n <k{\\Al' 2 i;\\ + U\\) <2*||A£/V||. 

This shows that the first GJN condition, (3.1), is satisfied for our triple. 
Next, we write 

(xApX^P, A P i>) = (A P x^,A pX ^) + Ri, (A.ll) 

where 

R 1 = (A vX <l>,\x,h v \<l>). (A.12) 
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Since x^P £ ^(Ap), an d A p is essentially selfadjoint on C^(M 3 ), there exists 
a sequence {y? n } C C^°(1R 3 ) s.t. y? n — > x - ^ an d A p y9 n ~" > ^-pX^P- Moreover A p 
leaves C^°(R 3 ) invariant, so we have 

n 

= \im{(A P xip,A p ip n ) + (xi/j, [A p , A p ]ip n )} , (A.13) 

n 

and we calculate (strongly on C£°(IR 3 )): i [A p , A p ] = H p + W - x ■ Vv - 2x 2 . 
Since x(H p + W — x ■ Vi>) is bounded, and x^P £ ^(a? 2 ) (because it is in 
£>(A P )), we get 

(A.13) = lim (A pX ^, A pVn ) + i? 2 , (A.14) 

n 

where we defined 

R 2 = x(H p + W - x ■ Vv - 2x 2 )x^) ■ (A. 15) 

Next, it is not difficult to see that if ip e C£° then x^ £ ^(A 2 ). Consequently, 
we can move A p in (A.14) to the left factor in the scalar product (recall that 
T>(A P ) D V(A P )), perform the limit, and move A p back to the right factor. 
We then obtain 

{A pX ip, A pX ip) = {A pX ip } A pX ip) + R 2 = (Apik, xA P x^) +R2-R1, (A. 16) 
where the bar denotes complex conjugate. Together with (A. 11) this gives 

{xA pX ^, A p i>) - (A P V, xApX^P) =Ri + R 2 ~ Ri- (A. 17) 

Let us first consider R 1 . We estimate 

< IIA^II I|[X, AplV'll < (IIA^II + \\[A p ,xM) lltxAM- 

It is clear that H-ApV'H < fc||Ap^ 2 ^||, and by the same argument as the one 
leading to (A. 10), that || < k\\ip\\ (use that x = F 2 > 0, as in the 

proof of Lemma A. 2), so that 

|^| < A;||A^|| ||[x,A p ]^|| < fc||AyV||(|MI + I|[X,^II), (A.18) 

where we used in the second step that [x, —A] is bounded. Next, Lemma A. 2 
(with X = H p , Z = x n , Y = A p , M = 1) shows that [x, x n ] is bounded, hence 
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\\[X,X 2 \n < En(k\\Xn4>\\ + \\x n [ X ,xM\) < k || || + En IWx, X n ]^ || . 

We write x n [x,x n ]ip = {x,x n }x n ip + [x n , [x,x n ]^. As above, || [x, x n ]x n ^\\ < 
fc||A]/ 2 ^||, and using x 2 — F, we write 

[x n , [x, x n ]] = -2[F, x n ] 2 + F[x n [F, x n ]] + [x n , [F, x n \\F. (A.19) 

As above, [F, x n ] is bounded, and a by now standard commutator expan- 
sion shows that so are the other two terms in (A.19). We conclude that 
||[x,a; 2 ]^|| < A;(||AyV|| + < fc||Ap /2 ^||. Combining this with (A.18) 

yields 

\Ri\ < fcllA^f. (A.20) 
Finally, let us obtain the same upper bound for i?2- Since x(H p + W — x- Vt> ) 
is bounded we only need to show that | (xi^, x 2 x^ ) ) I < ^||A^/ 2, ?/>|| 2 . Using that 
[x, x n ] is bounded, we arrive at 

i (x^x 2 x ^) i = \m\\ 2 < m\ 2 + k J2 \M\\ 2 < k\\K /2 ^\\ 2 , 

n n 

which shows that 

\R 2 \ < k\\Al /2 4>\\ 2 . (A.21) 

Combining (A.17), (A.20) and (A.21) shows that (xA p x, A P , C^(M 3 )) satis- 
fies the second GJN condition, (3.2). ■ 

A. 3 Proof of Proposition 3.4 

We give the following lemma without a proof, which is not difficult to find, 
e.g. by using the results of Appendix A.l. 

Lemma A. 4 Letv G C p ~ 1 (W i ) be s.t. x a d a v is bounded, for any multi-index 
\a\ < p — 1. Let H = —A + v, which is essentially self adjoint on C^°(lR <i ). 
Given a function /i G C£°(R) ; define fi(H) = f fi(s)e tsH . Then we have 

\\{x)* n n(H){x)± n \\<K(p,n), n = 0,l,...,p, (A.22) 

where K(p, ji) is some finite constant. 
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Lemma A. 5 The regularized G a> j = (pj d + fi)G a (pj d + /i) satisfies the same 
bounds (2.22) as G a . Moreover, adj^ (G a ,j) is bounded, n — 0, 1,2,3. 

Proof. The first assertion follows easily from the fact that (x) m pj d (x) n is 
bounded, for all m,n, and from Lemma A. 4 (use (x) n fi = (x) n fi(x)~ n (x) n ). 

In order to show boundedness of the multi-commutators, we treat a typ- 
ical term appearing in ad^ (Gr ajt /): 

xA p xG a ,jX A pXX A pX 
= xA p {x)~ l (x)x(x)' 1 (x)G aJ (x) 2 (x)~ 2 xA P x(x) (x)' 1 A P x- 

Since xA p {x)^ 1 and (x)^ 1 A P x are bounded, we see from Lemmas A. 4 and the 
bound (2.22) (forG aj j), that the r.h.s. is bounded, provided || {x)~ 2 xA P x(x} || < 
oo. To obtain the latter bound, it is enough (due to due to Lemma A. 4) to 
show || (x)~ 2 xA p (x) || < oo, which in turn is proved by writing 

{x)~ 2 X A p {x) = % -(x)- 2 xY,( X n(x)dn + 1/2), 

n 

and proceeding as in the proof of Lemma (A. 4), by commuting x n (x) through 
X to the left. ■ 

Proof of Proposition 3.4- The operator x{H p + W)x is bounded, hence 
relatively A p -bounded. We will use below the fact that [x, A p ] is bounded, 
which follows from Lemma A.2, with X = H p , Z = A p , Y = A p , M = 1. We 
have, in the strong sense on V(k p ), 

&d (2) ApX (H p ) = xlx 2 H p , A p ]x + {xWx, xA P xl 

where the commutator in the first term is bounded, and the second term 
equals x(Wx 2 A p — A p x 2 W)x, which is easily seen to be bounded, too (use 
Lemma A. 4 together with the fact that xApix)' 1 is bounded, and so is 
W(x)). Next, we show that 

*& { X X {H P ) = [x[x 2 H p , A P ] X , X A pX ] + [\XWX, xA pX ], xA pX ] (A.23) 

is bounded. The first term on the r.h.s. is the sum of two bounded operators 
plus x[ [x 2 H P) Ap], xA P x\Xi which can be written (by commuting x through 
A p ) as a bounded operator plus 

x[[x 2 H p ,A p },A pX 2 ]x- (A.24) 
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Setting xi = xH p , we expand \xi, A p \ = x[Hp + f Xi( s ) Jq e^ s ^ Sl ^ Hp We iSlHp , 
and (A. 24) splits into two terms, the first one, x[x[Hp, Aplx 4 , is bounded. 
To see boundedness of the second term, write it as 

/ L x{{ ei(S ~ Sl)HpWetSlHp i A P^ 2 + A P ei{S ~ S ^ 

and use that (x)W is bounded. The second term in (A. 23) can be written as 
a bounded operator plus x\Wx 2 A p — A P x 2 W, x 2 A p ]x- The latter term equals 
again a bounded operator plus 2Re[Wx 2 A p , x 2 A p ], which is easily seen to be 
bounded, by again noticing that (x)W, and derivatives of W multiplied by 
(a;) 2 are bounded. 

We have thus shown that the multi-commutators appearing in (3.42) are 
bounded (hence A p -bounded in the sense of Kato). Clearly, N is A/ bounded, 
and Lemma A. 5, together with the fact that ^p{{—id u ) k Tf}(g a ))) is relatively 
Af ^-bounded, shows that /„ is Aj-bounded in the sense of Kato. Conse- 
quently, condition (3.1) is satisfied for X = C n . 

Next, we verify that (3.2) is satisfied. Let us start with the commutator 
of Ci with A. We need to show relative boundedness of [x(H + W)x, A p ] and 
[ii,A] (relative to A p and Aj, respectively, in the sense of quadratic forms). 
Noticing that A p = H p — v + x 2 , we write [x(H + W)x, A p ] as a sum of a 
bounded operator plus 

[x 2 H p ,x 2 ] + [ X W X ,x 2 }. (A.25) 

Now setting xi — X 2 H P) the first term in (A.25) equals J2 n (x n [xi, x n ] + 
[Xi,Xn]x n ), so for any ip E C£°, 

\(^,iXi,x 2 ^)<kJ2\M\\ U\\<k(^,A p ^). (A.26) 

n 

Next, 

\(i>,\xW X ,x 2 }i>)\ (A.27) 
< 2| (^, X W[x,x 2 }^) | <2^|(^, ( X Wx n [x,x n } + X W[x,x n }xM |. 

n 

Commuting x n in the first term in the sum through x to the left, one sees that 
| {i/j,xWx n [x,x n ]i/j) | < fcdl^H 2 + ||x n ^|| \\ip\\), which is bounded from above 
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by k (ip , Apip) (proceed as in (A. 26)). The second term in the sum in (A. 27) 
is estimated in the same way. This shows that (A. 25) is A p -form-bounded. 

Next, in order to show the relative bound on [I n , A], it is enough to show 
that [ad^ x (G Q ), A p ] is relatively Ap-form-bounded, and that 

MHa.) n T fi ( ga ),AA 

is relatively A/-form-bounded. The former bound is easily obtained from 
(2.22), and the latter has been treated in subsection 3.3.1. This shows that 
I n are relatively A-form-bounded, hence also completing the proof that C\ 
satisfies condition (3.2). 

Next, we consider the commutator of C2 with A. The only thing to check 
is that [ad x A pX (Hp) , Ap] is A p -form-bounded. This commutator can be writ- 
ten as a bounded operator plus [ad^ A (H p ), x 2 }, hence it suffices to show 

that x n [ad^4 (H p ),x n ] is A p -form-bounded (n = 1,2,3). One shows that 

(2) 

[aa^ A (H p ), x n \ is bounded, by simple estimates as above. Relative bound- 

edness of x n [ad xApX (H p ) , x n ] then follows easily (proceeding as in (A. 26)). 
Consequently, (3.2) is satisfied for C 2 . 

We now consider the commutator of C3 with A, and it is enough to show 
that [a,d xApX (H p ) , x 2 ] is relatively A p -form-bounded. We write this commu- 
tator as 2Re [ad xApX (H p )xA P x, x 2 }. Now we have 

[ad%X x ( H p)xApX,x n ] = [^ ApX (H p ),x n ]xA p x + ad { x ] ApX (H p )[xA p x, x n ], 

and it is clear that [ad xApX (H p )xA P x, x n ] is bounded. Consequently, 

>, l^ ApX (H p ),x 2 ^)\ < km 2 + ty,*V» < a:||a;/^h 2 . 

Hence C 3 satisfies (3.2) and the proof of Proposition 3.4 is complete. ■ 
A. 4 Proof of Proposition 3.2 

We consider first the case when (2.20) holds. From UIU = we have 
UIR e IU = UIR 2 IU. We recall that Pq is the projection onto the kernel 
of L p and P = 11 — P is given by 

Po= ^2 Pm®Pn+Pd®Pc+Pc®Pd+Pc®Pc, (A.28) 

m,n(EA4 
E(m)y£E(n) 
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where p d and p c are the projections onto the discrete and continuous sub- 
spaces corresponding to H p . One can see that 

eUIR 2 £ P m^O, eUIR 2 (p m ®p n )m^0, 

E(m)=E(n) 

as e -> 0, that n/.R 2 (p c <g)p c )/n = 0, and that UIR 2 (p c <S>p d )IU = IlIR 2 (p d ® 
p c )IIl. From formula (2.45) for the interaction, we obtain the bound 

UIR 2 e m > TlIR 2 (p c ®p d )m 

= ^2 Z\2^ Pm ®Pn ® Pn){G a <g> tp <g> a(Tp(g a ))} 

E(m)=E(n)=E(m / ) 



We write a(Tp(g a )) = J RxS2 Tp(g a )(u,'E)a(u,'E) and use the pull through 
formula a(u, S)L = (L + u)a(u, S) and obtain 



-2 x - x _ rE(m) r 2 

n/^/n > ^ ^/ du dY, _ 0u _ - g a {-u,E)g^(-u,E) 

m,n,m'eM a,a> J 00 J S ' 2 

E(m)=E(n)=E(m') 



x ( p m G a — P ° 2 G a , Pm , )(S)P n (S) P n , (A.29) 

V {Hp - E(m) + u) 2 + e 2 J 

where we recall that E{m) is the eigenvalue of H p corresponding to the mode 
m. We have dropped the integration over the values u > E(m) because 
e((H p - E(m) + u) 2 + e 2 ) -1 -> 5(H P - E(m) + u) as e -> 0, hence it = 
—.Hp + E(m) < E(m). Recalling the definition of F, see before 2.19, and 
making the change of variable u \— > — w in the integral, we arrive at 



du / dE-^ (A.30) 

m,n,m' fcyw v, y 

E{m) = E(n) = E(m') 

x U m F(«,E) — — ^^F(u,E)*p m , ) ® Pn ®P^. 
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The projection p(E) onto the eigenspace corresponding to an eigenvalue E 
of H p is given by J2 meM . E( m )=EP™ and we use 

E = E E E E 

m,n,m'eM E£<T V (H P ) rneM neM m' eM 

E(m)=E(n)=E(m') E(m)=E E(n)=E E(m')=E 

to arrive at 




(A.31) 



yp(E)F(u, Z) (Hp _ E Z u)2 + e2 F(u, X)*p(E)) ® p(E) ® P n . 

The desired bound (3.27) now follows from (2.20) and (2.21). 

The case when (2.23) holds and 7 is given by (2.24) is done similarly. ■ 
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